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Abstract 

Numerical results for the distributions of light transmitted through metallic planar lenses com- 
posed of symmetric nanogroove arrays on the surfaces of a gold film are presented and explained. 
Both the near- and far-field distributions of the intensity of light transmitted are calculated by 
using a Green's function formalism. Results for an optimal transverse focus based on a quadratic 
variation of groove width are obtained. Meanwhile, a significant dependence of the focal length on 
the wavelength of light incident from the air side through the gold film into a dielectric substrate 
is found for this detector configuration. 



I. INTRODUCTION 



There is considerable interest at the present time in shaping the spatial dependence 
of the intensity of light transmitted through a metal film pierced by a one-dimensional 
array of slits with subwavelength widths, as an alternative to refractive lensing. Sun and 
Kim 1 have studied numerically the transmission of light through a finite periodic array 
of subwavelength slits piercing a free-standing metal film with a convex dependence of its 
thickness. The incident beam is focused by this metallic lens. Shi et al. ,? studied numerically 
the focusing of light transmitted through a metallic film of constant thickness pierced by 
an array of equally spaced nanoslits of varying widths. It was argued that the focusing 
action is based on propagation across the opaque metal film of surface plasmon polaritons 
supported by the slits (metal-dielectric-metal structures) of varying widths, thus ensuring an 
enhanced transmission and a phase change along the film surface that lead to the constructive 
interference. The focusing of light by such a structure deposited on a dielectric substrate 
was demonstrated experimentally by Verslegers et alM The experimental results were in 
an excellent agreement with FDTD simulations^. However, the focus distance was much 
shorter than predicted by the model proposed previously. 

In this paper we demonstrate the focusing effect by calculating numerically the spatial 
distribution of the intensity of light transmitted through a metal film sandwiched between 
a cladding and a dielectric substrate whose surfaces are modeled by two finite aligned and 
reversed arrays of nanogrooves of finite depth. We have chosen to define transmissivity as 
the squared modulus of the ratio of the transmitted to the incident H-field amplitude for 
p polarization. In this case there are no surface polaritons propagating through the film as 
there are no slits that completely pierce the metallic film. In the case of the extraordinary 
transmission of light through a metal film with a periodic nanoslit array, the surface polari- 
tons supported by the slits play a minor role in the phenomenon, and what is more it is 
not necessary to have slits that completely pierce the metallic film to achieve the enhanced 
transmission 4: . When a periodic nanogroove array is illuminated by p— polarized light, whose 
magnetic vector is parallel to the generators of the array, surface plasmon polaritons (SPPs) 
associated with the film-cladding and the film-substrate interfaces are excited^, and are 
diffracted by the structure into transmitted volume waves in the substrate. The periodic- 
ity of the array enhances the excitation of surface plasmon polaritons associated with the 
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cladding-film interface and the conversion of surface plasmon polaritons associated with the 
substrate-film interface into volume waves in a range of SPP frequencies whose wavenumbers 
are in the vicinity of the boundary of the second Brillouin zone, i.e. when X sp ~ d, where d 
is the period of the array and X sp is the surface plasmon polariton wavelength. As a result 
the transmission through gold and silver films as a function of wavelength^ shows sharp 
peaks and dips in a range of SPP frequencies whose wavenumbers are near the boundary of 
the second Brillouin zone at both interfaces of the film. 

It is also not necessary to have slits that completely pierce the metallic film to achieve 
the focusing of light transmitted through it when the array is finite. The conversion of 
either surface polaritons supported by the slits or surface polaritons supported by the film 
itself into the volume wave in the substrate is a diffraction process, and since the system 
is two-dimensional it produces volume cylindrical waves propagating away from the surface 
into the substrate. A finite number of secondary sources of cylindrical waves separated 
by subwavelength distances can produce the field focusing effect. The additional phases of 
each of the secondary sources changes the focal depth and width, but affect only slightly 
the focus distance. The effect can be modeled as follows. The periodic array of slits or 
nanogrooves ensures the transmission through the otherwise nontransparent film, while the 
now transparent aperture produces the focusing, in exactly the same manner as a refractive 
nanolens in a diffractive regime. The focus distance is then determined primarily by the 
aperture size, i.e. by the size of the array, and the wavelength of light in the substrate. The 
intensity distribution in the far field is then determined by the Fresnel diffraction by the 
aperture. 

In this paper we discuss focused patterns of nanogroove arrays and do not elaborate on 
differences between nanoslit and nanogroove arrays, since they are physically equivalent in 
producing the focusing effect except for some quantitative difference in the intensity of the 
transmitted field. A simulation study was undertaken to search for the best way to improve 
the quality of focus by changing the groove profile and groove width variation. Planar 
nanolenses will have numerous applications in polarimetric imaging devices, solar cells, light 
emitting diodes, and nanophotonics systems. Some preliminary results of this work was 
reported earlier.^ 

This paper is organized as follows. In Section II, we employ a previously developed 
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model and formalism,^ modified to include an arbitrary sequence of groove width variation. 
Based on this formalism, numerical results are presented for the comparison of the focused 
patterns of field intensity produced by different sequences of groove width variation and 
different groove shapes. The conclusions drawn from these results are briefly summarized in 
Section III. 

II. MODEL AND NUMERICAL RESULTS 

In this paper, we consider the model shown in Fig.[T]for a nanogroove array, in which the 
top [£i(x) — £(x) at z — 0] and bottom [^(x) = — £(x) at z = —L] surface profile functions^ 
employed for modeling a nanogroove array patterned in a thin gold film are chosen to be 

M 

£0) = E ex p 

j=-M 

where L is the thickness of the unpatterned gold film, j is the index for labeling grooves, 
2M + 1 is the total number of grooves in the array, t < L/2 is the groove depth, d is the 
period of the groove array, and {bj} is an arithmetic sequence representing a specific pattern 
of groove-width variation. The surface profile functions ^(x) in Eq. Q are symmetric with 
respect to the middle groove centered at x = 0, and s = 2 (or s = 4) corresponds to a 
Gaussian (or a quartic) functional form for a groove, respectively. The arithmetic sequence 
{bj} in Eq. (jl]) is assumed to be 

b J =a + P\ 3 \+ 1 \ 3 \\ for|j|<M, (2) 

where a represents the width of the central groove, and (3 = (or 7 = 0) corresponds to a 
quadratic (or a linear) groove-width variation, separately. The spatial distributions of the 
electromagnetic fields on the air side (n = n a ), inside the patterned gold film (n = n m ), and 
on the side of the dielectric substrate (n = n s ) can be calculated by using a Green's function 
formalism. SE1 

In our numerical calculations, whose results are presented below, we assume that the 
metal film is illuminated by a normally incident p— polarized plane wave, whose magnetic- 
field component H y (x, z) has a unit amplitude, or by a normally incident s— polarized plane 
wave, whose electric-field component E y (x, z) has a unit amplitude. The parameters defining 



_ ( x - jd \ 



4 



the film have the values L = 0.4/zm, t = 0.196 /im, d = 0.25 /im, M = 6, and a = 40nm 
for the width of the central groove. The frequency-dependent complex refractive index n m 
for the gold film is obtained by interpolation from the data in the paper by Johnson and 
Christy^. The values of the other parameters, such as n a , s , 0, 7, s, and the incident light 
wavelength An, used in our numerical calculations will be given in the figure captions. 



A. Focusing by a finite aperture 



In Fig.[2]we present color level plots of the intensity of the field of light \H y (x, z)\ 2 (left) 
of p and \E y (x, z)\ 2 (right) of s polarization transmitted through the gold film from the air 
cladding into the dielectric substrate. The surface profile functions are arrays of periodic 
grooves of the same half widths (/3 = 7 = 0). 

As can be seen from plot presented in Fig. [2] the interference pattern of the transmitted 
light of both polarizations exhibits a focal spot with the same focal distance but the intensity 
of the s— polarized field is seven orders of magnitude weaker. The additional interference 
maxima are of the same strength as the primary focal spot. Thus, in spite of the fact that 
the groove arrays considered in our calculations are finite (the aperture size is 1.5 /im), and 
the period is significantly smaller then the wavelength of light in air and in the substrate, 
the enhanced transmission in p— polarization is obvious. The intensity distribution of the 
primary focal spot is described by 



\H y (x, z) 



S 



c 2n s (x + V/2) 
sj\ \z + L\ 

2n s (x + V/2) 
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(3) 



where S(x) and C(x) are the Fresnel integrals, and V is the array length. The focal length 
(the position of the maximum of the intensity) is zq ~ n s D 2 ll/(28An), and is inversely pro- 
portional to the wavelength of light in the substrate Xq/u s . From the calculated \H y (x, z)\ 2 
(p-polarization) and \E y {x, z)\ 2 (s-polarization) with groove width variation (not shown 
here), we find that while the constructive interference is completely destroyed by introduc- 
ing the quadratic groove width variation in the case of s— polarized light, the quality of 
the primary focal spot is considerably improved and the subsidiary maxima of the inter- 
ference pattern are suppressed in the case of p— polarized light. Note that the focal length 



decreases. From the calculated total (integrated over the transmission angles from —tt/2 
to 7r/2) transmission coefficient as a function of the wavelength of the p— polarized incident 
light for different combinations of the cladding and substrate materials (not shown here), we 
clearly see the enhanced transmission in comparison with n a = n s = 1, although the array 
is limited to a few wavelengths. In addition, the introduction of a weak quartic change in 
the width of the grooves leads to a slight shift in the position and the broadening of the 
enhanced transmission peak, but does not destroy it. 

B. Wavelength effect on the focal length 

Figure [3] presents a comparison of the spatial distributions of \H y (x, z)\ 2 when 
p— polarized light is incident normally from the upper air side (corresponding to a detector 
configuration) at Ao = 0.63 /im (left panel) and Ao = 1 /im (right panel). Here, a quartic 
functional form (s = 4) is assumed for the symmetric nanogroove array on a gold film, which 
has a quadratic groove-width variation ((3 = 0). It is clear from Fig. [3] that the focal spot 
shifts upward from z = — 7/xm to z = —5 /im when Ao is increased from 0.63 fim to 1 /im. 
The positions of the maxima of the intensity are in agreement with the inverse dependence of 
the focal length on the wavelength. Meanwhile, the longitudinal [z direction) size of the fo- 
cal spot shrinks, although its transverse [x direction) size remains constant. This provides a 
possibility for multi-color detection if several active detection layers with specific absorption 
wavelengths are embedded at different depths underneath the nanogroove array on top of a 
quantum- well photodetector.E' The magnification also decreases, which is obviously related 
to the fact that the overall transmission decreases as the wavelength increases. 

In Fig.[4]we present the color level plot of the spatial distribution of the field intensity in 
the case when the dielectric substrate is adjusted to produce the strongest transmission at 
the given wavelength. In this case, the wavenumbers of SPP polaritons associated with the 
film-substrate interface are the same, so that at both wavelengths the effective aperture for 
SPP is the same. This is the reason for the same focal distance. It is clearly seen that the 
magnification is greatly increased for both chosen wavelengths. 
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C. Difference in detector and emitter configurations 

For device-configuration comparison, we show in Fig. [5] the results for \H y (x, z)\ 2 when 
the roles of air and dielectric substrate are interchanged, so that p— polarized light is incident 
normally from the upper cladding side with n a = 1.46 and is transmitted into the air n s — 1 
(corresponding to an emitter configuration) at Ao = 0.63 jum (left panel) and Ao = 1 //m 
(right panel), respectively. Here, a quartic functional form for a groove, as well as a quadratic 
groove-width variation, are also assumed for modeling the nanopatterned symmetric groove 
array. In this case, however, we find no shift of the focal spot in the z direction with a 
change of Ao- In this configuration the wavelength in the substrate (air) is considerably 
larger than in the case described in the preceding section. This leads to shorter focal length 
(1/1.46 times). Since the decay length ( of the groove near field in the dielectric substrate 
side scales like ( ~ Ao/n s ,^ the inverse dependence of the focal length on the wavelength 
will be completely masked by the increase of £ with Ao- At the same time, the transmitted 
field is reduced, as is seen by comparing Fig.|5]with Fig. [3] In addition, the transverse size 
of the focal spot is enlarged on the air side as Ao is increased from 0.63 /im to 1 /im. From 
the results shown in Fig.[5j we expect that this emitter configuration can be employed for 
constructing a wavelength-insensitive planar lens with a very thin dielectric substrate that 
is transparent to light of wavelength Ao > 0.63 /im. 

We also study the effect of a dielectric substrate on focusing action from the calculated 
\H y {x, z)\ 2 for p— polarized light incident normally from the upper air side at Ao = 1 /xm 
(not shown here), from which we find that the focusing power almost disappears when the 
dielectric substrate is replaced by air. Meanwhile, the non-focused bright spot is seen to 
move upward as n s is reduced from 1.46 to 1.0. On the other hand, the focal spot shifts 
downward as n s is increased from 1.0 to 2.0, in agreement with the diffractive regime of the 
lens, which is accompanied by an expansion of the longitudinal size of the focus spot. 

D. Effects of groove shape and groove- width variation 

To study the effect of groove shape, as well as the effect of groove-width variation, we com- 
bine four calculated spatial distributions for \H y (x, z)\ 2 in Fig. [6] obtained with p— polarized 
light of wavelength Aq = 0.8 /im incident normally in an emitter configuration. A quartic 
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functional form is assumed for the groove array (two upper panels) with a quadratic (left) 
and a linear (right) groove-width variation, to investigate the effect of groove-width variation. 
We also show the result for \H y (x, z)\ 2 (lower-right panel) with a constant groove width. In 
addition, the calculated \H y (x, z)\ 2 (lower-left panel) for the case with a Gaussian functional 
form (s = 2) is included in Fig.|6]to demonstrate the groove-shape effect. By comparing the 
two upper panels of Fig.[6j we find that the quadratic variation of groove width leads to an 
enhanced focusing power with a smaller spot size in both directions. For a constant groove 
width, the focusing action is completely lost. As the corners of a groove are rounded by a 
Gaussian functional form, in comparison with a quartic one, the focusing power is partially 
suppressed, but the focal spot does not move at all. 



III. CONCLUSIONS AND REMARKS 



In conclusion, when light is incident from the upper air side, we have demonstrated that 
the focal length of a planar metallic lens based on a variable nanogroove array deposited on 
a dielectric substrate can be controlled by varying the wavelength of the incident light. How- 
ever, this wavelength-tunability of the focal length for a metallic planar lens is completely 
absent when the light is incident from the substrate side. These numerical results can be 
applied to the design of a multi-color photodetector in which a number of active detection 
layers with specific absorption wavelengths are embedded at different depths underneath 
the top groove array. Moreover, the enhanced focusing power of a metallic planar lens with 
a quadratic groove-width variation is observed in comparison with that of a linear one. A 
range for the refractive index of a dielectric substrate is found as an imposed restriction 
for the tunable focal length of a metallic planar lens. The sharpness of a groove corner is 
shown to play an important role in the focusing power of a metallic planar lens through the 
comparison of a groove in a quartic functional form with that in a Gaussian one. 

When a p-polarized incident light illuminates a periodic nanogroove array, SPPs, which 
are associated with the film-cladding and film-substrate interfaces, can be excited. These 
excited SPP waves will be diffracted by the array into volume cylindrical waves after their 
transmission as long as the SPP wavenumbers fall into the vicinity of the boundary of the 
second Brillouin zone. Diffraction induces a coupling between SPP modes with different 
reciprocal lattice vectors and produces a gap at either the center or the boundary of the first 
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Brillouin zone. These generated cylindrical waves, which are separated by subwavelength 
distances, lead to the field focusing effect. We note that the periodic array of nanogrooves 
is only responsible for the transmission through the otherwise nontransparent film. It is the 
transparent aperture that produces the focusing. As a result, the focus distance is decided 
by both the size of the array and the wavelength of light in the substrate, which allows for 
a tunable focusing. 

When the total number, 2M + 1, of grooves in a metallic planar lens is increased, we 
expect to see a significant reduction in the transverse size of a focus spot, which facilitates 
an even smaller pixel size for a photodetector focal-plane-array to improve its detectivity, 
as well as a suppression of cross-talk between different pixels.^ 
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FIG. 1: (Color online) The structure studied, where n is the index of refraction, the spatial 

directions, £(x) is the surface profile function, and L is the thickness of the unpatterned gold film. 
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FIG. 2: (Color online) Color level plot of the spatial distribution of the intensity of the p— (left 
panel) and s— polarized (right panel) light of Ao = 0.63 fim transmitted through the gold film. A 
quartic functional form of the same width (/3 = 7 = 0) is used for the grooves constituting the 
surface profiles. The medium of incidence is air n a = 1, while the substrate is a dielectric with the 
refractive index n s = 1.46. 
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FIG. 3: (Color online) Contour plot of the spatial distribution of \H y (x, z)\ 2 for two incident-light 
wavelengths, Ao = 0.63 /xm (left panel) and Ao = 1 (right panel). Here, a quartic functional 
form is used and the light is incident from the upper air side (detector configuration). The other 
parameters in the calculations are n s = 1.46, as well as f3 = and 7 = 35/36 nm for the quadratic 
(or parabolic) groove-width variation. 
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FIG. 4: (Color online) Contour plot of the spatial distribution of \H y (x, z)\ 2 for two incident- 
light wavelengths and different substrates Ao = 0.63 //m , n s = 1.56 (left panel) and Ao = 1 /mi, 
n s = 2.61 (right panel). Here, a quartic functional form is used and the light is incident from the 
upper air side (detector configuration). The other parameters in the calculations are (3 = and 
7 = 35/36 nm for the quadratic groove-width variation. 
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FIG. 5: (Color online) Contour plots of \H y (x, z)\ 2 at Ao = 0.63 /im (left panel) and Ao = 1 //m 
(right panel). Here, a quartic functional form is assumed and the light is incident from the upper 
substrate side (emitter configuration). The other parameters are n a = 1.46, n s = 1, as well as 
/3 = and 7 = 35/36 nm for the quadratic groove-width variation. 
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FIG. 6: (Color online) Contour plots of \H y (x, z)\ 2 with the quadratic groove-width variation 
(upper-left panel) and the linear groove- width variation (upper-right panel). Here, a quartic func- 
tional form is assumed for the upper two panels in this figure. We take (3 = and 7 = 35/36 nm 
for the quadratic groove- width variation and f3 = 35/6 nm and 7 = for the linear groove- width 
variation, respectively. We also display here the contour plots for the scaled \H y (x, z)\ 2 with the 
quadratic groove-width variation (lower-left panel) and with a constant groove width (lower-right 
panel). Here, a Gaussian functional form is used for the lower- left panel, while a quartic functional 
form is assumed for the lower-right panel. In addition, we take (3 = and 7 = 35/36 nm for the 
quadratic groove-width variation and (3 = 7 = for the constant groove width, separately. The 
light is incident from the upper substrate side at Ao = 0.8 /xm, and the refractive index of the 
substrate is n s = 1.46, where n a = 1 and the black curves close to z = indicate the surface 
profiles of the patterned gold film. 
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